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Abstract 

It was a difficult problem to determine the Gaussian fixed line 
from the numerical data, because close to the Berezinskii-Kosterlitz- 
Thouless multicritical point the divergence of the correlation length 
becomes very slow. Considering the renormalization group behavior, 
we find an efficient method to determine the Gaussian fixed line. 

PACS numbers: 75.40. Cx, 05.70.Fh, ll.10.Hi, 75.10.Jm 



1 



The two dimensional (2D) sine-Gordon model, which is defined with the 
Lagrangian 



plays an important role in 2D classical and ID quantum systems, such as 
the 2D XY model, 2D helium, ID quantum spin models and ID electron 
models. It is also related with the U(l) x Z 2 x Z 2 conformal field theory or 
the Tomonaga-Luttinger liquid. The sine-Gordon model has the Berezinskii- 
Kosterlitz-Thouless (BKT) 0, [| || phase boundaries and the Gaussian fixed 
line. 

However, it was a difficult problem to determine the phase boundary and 
the universality class of the sine-Gordon model from the numerical calcula- 
tion. Finite-size scaling technique, which is an efficient method to analyze 
the ordinary second order transition, leads to false conclusions for the BKT 
transition [|], [5|, |6], [7|, since the divergence of the correlation length is very 
slow £ ~ exp(const./ \/T — T c ) and there are logarithmic corrections caused 
by the marginally irrelevant field. One of us || has proposed a remedy for 
these problems, the level spectroscopy, based on the renormalization analysis 
and the SU(2)/Z2 symmetry of the BKT transition. 

About the Gaussian fixed line, in the neighborhood of the BKT multi- 
critical point, since the divergence of the correlation length is very slow, it is 
difficult to determine the critical line with the finite-size scaling technique. In 
addition, from the simple finite-size scaling analysis of the Gaussian model, 
two (pseudo) fixed points or no fixed point are obtained, which causes a 
finger- like massless region around the true Gaussian fixed line |9|, 0, [□], |]. 
In this letter we propose an improved method to obtain the Gaussian fixed 
line. 

The renormalization group equations for ([I]) are || 



(with the identification = + 27r/v2) and the dual field 9 




(1) 



40 = -d y (iK6), dycf) = d x {iK6), 



(2) 



dl 
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For the finite system, I is related to L by e l = L. Close to the BKT multi- 
critical point K « 4, \y\\ <C 1, the divergence of the correlation length close 
to the Gaussian line becomes very slow, as mentioned above. 

On the Gaussian fixed line (yi = 0), the operators O n ^ m = exp(m\/20 + 
im\/29) has a critical dimension x n ^ m and a spin l n ^ m given by |L2| 



x n ,m = ^ (n 2 K + , l njTn = nm. (4) 

When we denote the transfer matrix of a strip of width L with periodic 
boundary condition by exp(— H), then the energy gaps AE n)in are related to 
the scaling dimension as [|T3[ 

AE n , m (L) = (5) 

where v is the "velocity of light" . Therefore, we obtain 

AE 0t2 (L)/AE 0A (L)=4. (6) 

On the other hand, the off-critical behavior y\ ^ 0, K < 4 is considered 
as follows. In this case eq. (|f) is infrared unstable and the correlation length 
£ becomes finite. In the L>( limit, it is considered that the excitation of 
m = 2 is a scattering state of two excitations of m = 1, which means 

AE 0>2 (L)/AE 0jl (L) = 2. (7) 

Therefore, we expect that the ratio AEo^{L) / AEqi(L) becomes maximum 
at yi — and decreases as increasing \yi\ (see Figure 1). 

Consider the other type ratio AEq^{^L) / AE ^(L) . It takes 2 at the 
Gaussian line, and in the massive region it also takes 2, therefore it is expected 

A£ , 2 (2L)/A£ ,i(£) = 2 (8) 

independently of yi (see Figure 2). 

We shall explain these behaviors with the renormalization group argu- 



ment. The renormalized critical dimensions are |p~4" | 

-,2 



x , m = —K-\l) (9) 
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However, with only these relations, we cannot explain the behavior AEq^{L) / AEq^{L) 
close to the Gaussian line, since the ratio of xo, m is independent of y±. Con- 
sidering the renormalization constant L , which is included in the definition 
I as I = log(L/L ), we can understand these features (6), (7), (8). In order to 
satisfy the relation @, the renormalization constant for the AEq^ should be 
2L , where L is the renormalization constant for AE 01 . Then, we obtain 



AE , 2 (L) 



K-\l - log 2) 
K-\l) 



K-\l) 



K-\l)-\og2 



- l=logL/Lo 

dK-\l) 



dl 



- Z=logL/Lo 



log 2 



l-jr(Q-2-yf(J) 



(10) 



l=logL/L 



by using eqs.(^). Consequently we can determine the Gaussian fixed line 
from the extremum point of the AE 0)2 (L) j AE 0i i(L) . 

As a physical model, we study the S=l bond-alternating XXZ model [ITS 



H 



L 

B 



i 



5(-iy)(S*S* +1 + S]S] +1 + ASjS, 



3+1' 



(11) 



In this case AE ^ m (L) corresponds to AE(S^ = m,L), where Sj, = J2Sj. 
In Figure 1, we show the ratio AE(2, L)/AE(1, L). And in Figure 2, we 
show the ratio AE(2, 2L)/AE(1, L). Their behaviors are consistent with the 
renormalization discussion. The remaining correction 1/L 2 can be explained 
by the irrelevant operator L_ 2 -^-2l with the critical dimension x = 4. 

Here we analyze the corrections caused by the irrelevant operator L_ 2 X_ 2 1- 
For simplicity we treat the case \yi\ <^ \K—A\. In this case y\ = yi(0)L 2 ~ K ^ 2 . 
And we assume the correction from the x = 4 term as (c x + c 2 (5 — 5 C ))L~ 2 
and the dependence of y\ on the parameter 5 as y±(0) = c 3 (5 — 5 C ). Then the 
extremum of the AEq 2 {.L)/AEq i{L) is situated at 



5-5 r 



C'2 



K-6 



K\og2cl 



(12) 



which rapidly converges to in the L — > oo limit. 

Finally we criticize the method proposed in |_ 
appears the term y?(l)/(K(l) — 4). However, close to the BKT multicritical 



In this method, there 
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point, the sign of K{1) — 4 changes, therefore it becomes difficult to determine 
the maximum. 
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Figure 1: Ratio of excitation energies AE(2, L)/AE(1, L) for L — 10 (x), 
L — 12 (□), L = 14 (+) and L — 16 (O). (a) A = where the critical value 
of 5 is 5 = 0.23. (b) A = 0.5 where the critical value of 5 is 5 = 0.25. 
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Figure 2: Ratio of excitation energies AE(2, 2L)/AE(1, L) for L — 4 (x), 
L = 6 (□) and L = 8 (+). O is th^j extraporated value, (a) A = 0. (b) 
A = 0.5. 



